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Abstract: A Nielsen-Olesen vortex usually sits in an environment that expels the flux that 
is confined to the vortex, so flux is not present both inside and outside. We construct vortices 
for which this is not true, where the flux carried by the vortex also permeates the ‘bulk’ far 
from the vortex. The idea is to mix the vortex’s internal gauge flux with an external flux 
using off-diagonal kinetic mixing. Such ‘dark’ vortices could play a phenomenological role 
in models with both cosmic strings and a dark gauge sector. When coupled to gravity they 
also provide explicit ultra-violet completions for codimension-two brane-localized flux, which 
arises in extra-dimensional models when the same flux that stabilizes extra-dimensional size 
is also localized on space-filling branes situated around the extra dimensions. We derive 
simple formulae for observables such as defect angle, tension, localized flux and on-vortex 
curvature when coupled to gravity, and show how all of these are insensitive to much of the 
microscopic details of the solutions, and are instead largely dictated by low-energy quantities. 
We derive the required effective description in terms of a world-sheet brane action, and 
derive the matching conditions for its couplings. We consider the case where the dimensions 
transverse to the bulk compactify, and determine how the on- and off-vortex curvatures and 
other bulk features depend on the vortex properties. We find that the brane-localized flux 
does not gravitate, but just renormalizes the tension in a magnetic-field independent way. 
The existence of an explicit UV completion puts the effective description of these models on 
a more precise footing, verifying that brane-localized flux can be consistent with sensible UV 
physics and resolving some apparent paradoxes that can arise with a naive (but commonly 
used) delta-function treatment of the brane’s localization within the bulk. 
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1. Introduction 

In this paper we study the gravitational response of vortices that carry localized amounts of 
external magnetic flux; called Dark Strings or Dark Vortices in the literature [1, 2]. The goal 
is to understand how their back-reaction influences the transverse geometry through which 
they move, and the geometry that is induced on their own world-sheet. We find the initially 
surprising result that the gravitational response of such an object is locally independent of 
the amount of flux it contains, and show how this can be simply understood. 
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Motivation 

Why study the gravitational response of Dark Vortices? 

Vortices are among the simplest stable solitons and arise in many theories with sponta¬ 
neously broken 17(1) gauge symmetries [3]. They can arise cosmologically as relics of epochs 
when the Universe passes through symmetry-breaking phase transitions. Such cosmic strings 
are widely studied [4] because, unlike other types of cosmic defects, they need not be poi¬ 
sonous for later cosmology since the resulting cosmic tangle tends not to come to dominate 
the energy density in a problematic way. 

In the simplest models a vortex defines a region outside of which the 17(1) symmetry 
breaks while inside it remains (relatively) unbroken, and as a result all magnetic 17(1) flux 
is confined to lie completely within the vortex interior. However in theories with more than 
one 17(1) factor more complicated patterns can also exist, for which magnetic fields outside 
the vortex can also acquire a localized intra-vortex component. Such vortices naturally arise 
in ‘Dark Photon’ models [5], for which the ordinary photon mixes kinetically [6] with a 
second, spontaneously broken, 17(1) gauge held (as have been widely studied as Dark Matter 
candidates [7]). Cosmic strings of this type could carry localized ordinary magnetic flux, even 
though the U(1) E m gauge group remains unbroken [1, 2], 

Of most interest are parameters where the vortex’s transverse thickness is much smaller 
than the sizes of interest for the geometry transverse to the source. In such situations only 
a few vortex properties are important, including the tension (energy per unit length) and 
the amount of flux localized on the vortex (or more generally brane-localized flux, or BLF 
for short). Indeed these two quantities (call them 7), and Cb) provide the coefficients of the 
leading terms in any derivative expansion of a vortex action (for which more explicit forms 
are also given below), 

Sb = —Tb J to + Cb J * A + • • ■ , (1-1) 

where w is the volume form of the codimension-two surface and *A is the Hodge dual of the 
U( 1) field strength, A MN = d M A N — d N A M whose flux is carried by the vortex. These are the 
leading terms inasmuch as all terms represented by the ellipses involve two or more deriva¬ 
tives. 1 In four dimensions both u and *A are 2-forms and so can be covariantly integrated 
over the 2-dimensional world-sheet of a cosmic string, while in D = d + 2 dimensions they are 
d forms that can be integrated over the d -dimensional world volume of a codimension-2 sur¬ 
face. 2 Previous workers have studied gravitational response in the absence of brane-localized 
flux [8], but our particular interest is on how (b competes with Tb to influence the geometry. 
Our analysis extends recent numerical studies [2] of how dark strings gravitate, including in 
particular an effective field theory analysis of the BLF term and its gravitational properties. 

Besides being of practical interest for Dark Photon models, part of our motivation for this 
study also comes from brane-world models within which the familiar particles of the Standard 

single-derivative term involving the world-sheet extrinsic curvature is also possible, but our focus here 
is on straight motionless vortices. 

2 That is, a brane with precisely two transverse off-brane dimensions. 
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model reside on a 3+1 dimensional brane or ‘vortex’ within a higher-dimensional space. 3 
Comparatively little is known about how higher-codimension branes situated within compact 
extra dimensions back-react gravitationally to influence their surrounding geometries, 4 and 
codimension-2 objects provide a simple nontrivial starting point for doing so. In particular, a 
key question in any such model is what stabilizes the size and shape of the transverse compact 
dimensions, and this is a question whose understanding can hinge on understanding how the 
geometry responds to the presence of the branes. Since long-range inter-brane forces vary only 
logarithmically in two transverse dimensions, they do not fall off with distance and so brane 
back-reaction and inter-brane forces are comparatively more important for codinrension-2 
objects than they are with more codimensions. 

Furthermore, several mechanisms are known for stabilizing extra dimensions, and the 
main ones involve balancing inter-brane gravitational forces against the cost of distorting 
extra-dimensional cycles wrapped by branes or threaded by topological fluxes [11, 12, 13, 14, 
15]. Since brane-localized flux is the leading way fluxes and uncharged branes directly couple 
to one another, it is crucial for understanding how flux-carrying vortices interact with one 
another and their transverse environment. Localized flux has recently also been recognized 
to play a role in the stability of compact geometries [16]. 

Finally, the fact that cosmic strings can have flat world-sheets for any value of their string 
tension [8] has been used to suggest [17, 18] they may contain the seeds of a mechanism for 
understanding the cosmological constant problem [19]. But a solution to the cosmological 
constant problem involves also understanding how the curvature of the world-sheet varies as 
its tension and other properties vary. This requires a critical study of how codimension-2 
objects back-react onto their own induced geometry, such as we give here. Although extra- 
dimensional branes are not in themselves expected to be sufficient to provide a solution 
(for instance, one must also deal with the higher-dimensional cosmological constant), the 
techniques developed here can also be applied to their supersymmetric alternatives [20], for 
which higher-derivative cosmological constants are forbidden by supersymmetry and whose 
ultimate prospects remain open at this point. We make this application in a companion paper 
[ 21 ]. 

Results 

Our study leads to the following result: brane-localized flux does not gravitate. It is most 
intuitively understood when it is the dual field F = * A that is held fixed when varying the 
metric, since in this case the BLF term S BLF = ( f F is metric-independent. We show how 
the same result can also be seen when A is fixed; and more precisely show that the Q, (or 
BLF) term of (1.1) induces a universal renormalization of the brane’s tension and the brane 
gravitational response is governed only by the total tension including this renormalization. 

4 Our restriction to codimension-2 branes makes d = 4 and D = 6 the most interesting case of this type [9]. 

4 By contrast, back-reaction is fairly well-explored for codimension-1 objects due to the extensive study of 
Randall-Sundrum models [10]. 
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This renormalization is universal in the sense that it does not depend on the size of any macro¬ 
scopic magnetic field in which the vortex may sit. (The central discussion, with equations, 
can be found between eqs. (3.46) and (3.50), below.) 

Of course the BLF term does contribute to the external Maxwell equations, generating 
a flux localized at the vortex position with size proportional to Among other things this 
ensures that a test charge that moves around the vortex acquires the Aharonov-Bohm phase 
implied by the localized flux. But its gravitational influence is precisely cancelled by the 
back-reaction of the Maxwell held, through the gravitational held set up by the localized flux 
to which the BLF term gives rise. Since an external macroscopic observer cannot resolve 
the energy of the vortex-localized BLF term from the energy of the localized magnetic held 
to which it gives rise, macroscopic external gravitational measurements only see their sum, 
which is zero. 

The presence of the localized energy in the induced magnetic held does change the total 
energy density of the vortex, however, which can be regarded as renormalizing the vortex 
tension. This renormalization is independent of the strength of any outside magnetic helds. 

This failure of the BLF term to gravitate has important implications for the curvature 
that is induced on the vortex world-sheet. To see why, consider the trace-reversed Einstein 
equations in D = d + 2 dimensions, which state 5 


RmN T K ( TmN d 9mn T 


= 0 . 


( 1 . 2 ) 


What is special about this equation is that the factor of 1 /d ensures that the on-brane stress- 
energy often drops out of the expression for the on-brane curvature, which is instead governed 
purely by the off-br&ne stress energy. Consequently it is of particular interest to know when 
T m n vanishes for some reasonable choice of brane lagrangian. 

T mn would vanish in particular when the brane action is dominated by its tension 

rji _ rji n S{V) ( 

1 flU — ^b9fj,u J (T-'jj 

y/92 

where 5(y) is some sort of regularized delta-like function with support only at the brane 
position. But the derivation of (1.3) from (1.1) is complicated by two issues: is there a 
dependence on the transverse metric hidden in the regularized 5(y) (which is designed, after 
all, to discriminate based on proper distance from the vortex); and (for flux-containing branes) 
what of the metrics appearing in the Hodge dual, *A, of the BLF term? 

The results found here imply these two issues are not obstructions to deriving (1.3) from 
(1.1). They do this in two ways. First they show how T mn can be derived without ad-hoc 
assumptions about the metric-dependence of 5(y). Second, they show that the apparent 

5 We use Weinberg’s curvature conventions [22], which differ from those of MTW [23] only by an overall 
sign in the definition of the Riemann tensor. Coordinates x M label all D dimensions while x ^ ( x m ) label the 
d-dimensional (2-dimensional) subspaces. 
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dependence of the BLF terms on the transverse metric components, g mn , is an illusion, 
because it is completely cancelled by a similar dependence in the gauge-held back-reaction. 

The remainder of this paper shows how this works in detail. We use three different 
techniques to do so. 

• The first works within a UV completion of the dark vortex, for which we explicitly solve 
all held equations for a system that allows Nielsen-Olesen type vortex solutions. In this 
construction the BLF term can arise if there is a kinetic mixing, eZ MN A MN , between 
the 17(1) gauge held, Z M , of the Nielsen-Olesen vortex, and the external gauge held, 
A M , whose hux is to be localized. In this case the mixing of the two gauge helds can 
be diagonalized explicitly, leading to the advertised cancellation of the BLF coupling as 
well as a renormalization of the Z M gauge coupling, e 2 —> e 2 = e 2 /(l — s 2 ). 

• Second, we compute the couplings T and £ of the effective action for the codinrension-2 
vortex in the limit where the length scales of the transverse geometry are much larger 
than the vortex size. This has the form of (1.1), with ^ oc £/e. We verify that it 
reproduces the physics of the full UV theory, including in particular the cancellation of 
BLF gravitational interaction and the renormalization of the brane tension quadratically 
in £. 

• Finally we compare both of these approaches to explicit numerical calculations of the 
metric-vortex profiles as functions of the various external parameters, to test the ro¬ 
bustness of our results. 

A road map 

The remainder of the paper is organized as follows. 

The next section, §2, describes the action and field equations for the microscopic (or UV) 
system of interest. §2.1 shows this consists of a ‘bulk’ sector (the metric plus a gauge field, 
A m ) coupled to a ‘vortex’ sector (a charged scalar, T, and a second gauge field, Z M ). The 
vortex sector is designed to support Nielsen-Olesen vortices and these provide the microscopic 
picture of how the codimension-2 objects arise. The symmetry ansatze used for these solutions 
are described in §2.2 and the order-of-magnitude scales given by the parameters of the system 
are summarized in §2.3. 

Solutions to the field equations describing a single isolated vortex are then described in 
detail in §3, including both analytic and numerical results for the field profiles. The logic of 
this section, starting in §3.1, is to integrate the field equations in the radial direction, starting 
from initial conditions at the centre of the vortex and working our way out. The goal is to 
compute the values of the fields and their first derivatives just outside the vortex. In general 
we find a three-parameter set of choices for initial conditions (modulo coordinate conditions), 
that can be taken to be the flux quantum, n, for the vortex together with two integration 
constants (Q and R ) that describe the size of the ambient external magnetic field and the 
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curvature of the on-brane directions. 6 The resulting formulae for the fields and derivatives 
external to the vortex provide the initial data for further integration into the bulk, and are 
efficiently captured through their implications for the asymptotic near-vortex form of the bulk 
solutions, described in §3.3. In §3.4 these expressions for the near-vortex fields and derivatives 
are also used to match with the effective vortex description of (1.1) to infer expressions for 
Tb and (b in terms of microscopic parameters. 

The point of view shifts in §4 from the perspective of a single vortex to the question of 
how the bulk responds once the two vortices at each end are specified.' This is done in two 
ways. One can either continue integrating the field equations radially away from the first 
source (with n, Q and R specified as initial data as before) and thereby learn the properties 
of the source at the other end of the transverse space (by studying the singularities of the 
geometry where it closes off and compactifies). Alternatively, one can take the properties of 
the two sources as given and instead infer the values of Q and R that are consistent with the 
source properties: the two flux quanta n+ and n_, and the overall quantum N for the total 
magnetic flux in the transverse dimensions. After §4.1 first provides a set of exact integral 
expressions for quantities like R in terms of other properties of the source and bulk solutions, 
§4.2 describes the exact solutions for the bulk that are maximally symmetric in the on-brane 
directions and interpolate between any pair of source vortices. 

Finally, §5 summarizes our results and describes several open directions. Some useful but 
subsidiary details of the calculations are given in several Appendices. 

2. The system of interest 

We start by outlining the action and field equations for the system of interest. Our system 
consists of an Einstein-Maxwell system (the ‘bulk’) coupled to a ‘vortex’ — or ‘brane’ - 
sector, consisting of a complex scalar coupled to a second U( 1) gauge field. For generality 
we imagine both of these systems live in D = d + 2 spacetime dimensions, though the most 
interesting cases of practical interest are the cosmic string [with (D , d) = (4,2)] and the 
brane-world picture [with ( D,d ) = (6,4)]. 

2.1 Action and field equations 

The action of interest is S = S B + S v with bulk action given by 

^2 5 A/JV Kmn + ^ A M n^ IN + A 

=• — J d^ +2 a: y/—g ^ L EH + L A + A^ (2-1) 

6 For a given vortex lagrangian the tension of the vortex is controlled in terms of n by parameters in the 
lagrangian. We can also take the tension to be a separate dial - independent of n — if we imagine having 
several vortex sectors with different coupling constants in each sector. 

' Using electrostatics in 3 spatial dimensions as an analogy, §3 does the analog of relating the coefficient of 
1/r in the electrostatic potential to the charge defined by the properties of the source. Then §4 asks what the 
equilibrium configuration is for a coilection of charges given the resuiting electrostatic potential. 



- 6 - 





where A MN = d M A N — d N A M is a D-dimensional gauge field strength, 1Z MN denotes the D- 
dimensional Ricci tensor and the last line defines the L* in terms of the corresponding item 
in the previous line. The vortex part of the action is similarly given by 


Sv = - J a d+ 2 xV=g 

=: — J d d+2 x \J—g (^L z + L m ; x + Lq, + Vhj 


- Z mn Z mn + - Z mn A mn + + A [ T* T - V 


,. 2 \ 2 ' 


( 2 . 2 ) 


where := d M 'P — ieZ M 'F, and the second line again defines the various Li. 

For later purposes it is useful to write \/2 ^ and adopt a unitary gauge for which 

the phase, O, is set to zero, though this gauge will prove to be singular at the origin of the 
vortex solutions we examine later. In this gauge the term Lq, in S v can be written 

L* = DmT*^ (d M i> d M *l; + e 2 i> 2 Z M Z M ^J (2.3) 

and the potential becomes 

V b = - A (^ 2 -v 2 ) 2 . (2.4) 

It is also useful to group the terms in the brane and bulk lagrangians together according to 
how many metric factors and derivatives appear, with 


^kin := ^ 9 MN d M 9*rip , 


Lpot A + V b 


and 


Lgge : — Lj\ + L z + L m i x 

1 2 , 1.2 mn I 


L gm :=-e z ^g MN Z M Z N . 


Jgm • 2 

For this system the field equations for the two Maxwell fields are 


1 


V=g 


d h 


-g(A MN + eZ* 


and 


1 


V=9 


d a, 


-g [Z MN +e A* 


= 0, 


= e 2 ^ 2 Z N . 


The scalar field equation in unitary gauge becomes 
1 


V^~9 


d M [y~9 9 MN d N i= e 2 il>Z M Z M + A ^{^ 2 - u 2 ) , 


while the Einstein equations can be written in their trace-reversed form 


(2.5) 


( 2 . 6 ) 

(2.7) 


( 2 . 8 ) 




= -k 2 X 


MN i 


(2.9) 


where X MN := T MN — (1/d) T g MN and the stress-energy tensor is 


Tmn — d M i> d N ip + e 2 ip 2 Z M Z N + A mp A n p + Z MP Z N } 
£ 

2 


( 2 . 10 ) 


+- ( A mp Z n p + z 


■ M pA N ^ gMN ^ L]. j n -|- L gm T Lpot T Lgge^j 
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2.2 Symmetry ansatze 


We seek vortex solutions for which the brane/vortex sector describes energy localized along 
a time-like d-dimensional subspace, with nontrivial profiles in the two transverse dimensions. 
Accordingly, our interest is in configurations that are maximally symmetric in the d dimen¬ 
sions (spanned by and axially symmetric in the 2 ‘transverse’ dimensions (spanned by 

y m )• 

We take the fields to depend only on the proper distance, p, from the points of axial 
symmetry, and assume the only nonzero components of the gauge field strengths lie in the 
transverse two directions: A mn and Z mn . We choose the metric to be of warped-product form 


ds 2 = g MN dx M dx N = g rnn dy m dy n + g, JiU dx^dx" , (2.11) 

with 

gmn = 9mn(y) and g^ = W 2 (y) g^x ), (2-12) 

where g is the maximally symmetric metric on d-dimensional de Sitter, Minkowski or anti-de 
Sitter space. The corresponding Ricci tensor is 1Z M n dx M dx N = 1Z^ W dx^dx 1 ' + 1Z mn dy m dy n , 
and is related to the Ricci curvatures, R^ and R mn , of the metrics g liv and g mn by 


= R^ + g mn 


(d -1 )d m wd n w + wv m v n w 


9n v i 


(2.13) 


and 

R-mn = Rmn + V, ra V n W , (2-14) 

where V is the 2D covariant derivative built from g rnn . We work with axially symmetric 2D 
metrics, for which we may make the coordinate choice 


gmn dy m dy n = A 2 (r) dr 2 + B 2 {r) d 6 2 = dp 2 + B 2 {p) d6 2 , (2.15) 


where the proper radial distance satisfies dp = A(r)dr. With these choices the field equation 
simplify to the following system of coupled nonlinear ordinary differential equations. 

Gauge fields 

The gauge field equations become 


fw d A' e 


= 0 , 


(2.16) 


1 — £ 2 fW d Z' e V _e 2 *l) 2 Z e 

BW d V B ) ~ B 2 ’ 


(2.17) 


where primes denote differentiation with respect to proper distance, p, and we define the 
mixed gauge field, 


A m ■ — A m T £ Zm - 


( 2 . 18 ) 



Notice that the off-diagonal contribution to L gge vanishes when this is expressed in terms 
of A m rather than A M , since 


Laae — La + L z T- A m ; x — La + L 


Z > 


(2.19) 


where 

L A :=\A mn A mn and L z := i (1 - e 2 )Z mn Z mn . (2.20) 

Notice also that (2.17) has the same form as it would have had in the absence of the A — Z 
mixing, (2.18), provided we make the replacement e 2 —> e 2 , with 



( 2 . 21 ) 


Clearly stability requires the gauge mixing parameter must satisfy e 2 < 1 and semi-classical 
methods require us to stay away from the upper limit. 

Scalar field 

The field equation for fi’(p) similarly simplifies to 


1 

BW d 



( 2 . 22 ) 


Einstein equations 

The nontrivial components of the matter stress-energy become 

T^ = -g^Q t ot , T p p = Z — X and T e e = -{Z + X), (2.23) 


where 


and 


Q 7/kin + A gm + Lpot + L, 


'— Tp 0 t Lgg e 


(2.24) 


(2.25) 


and Z — Tkm T gm . 

In later sections it is useful to split g = qi oc + g B , X = X\ oc + X B and Z = Z\ oc + Z B into 
vortex and bulk parts, which we do as follows: 


6b •— A + La , £?loc •— Akin T Lg m -f- Vh + L^ 

X B := A — L a , X\ oc := 14 — L z 

and Z B . — 0 , -^loc •— Akin Lgm — E . 


(2.26) 


The components of the trace-reversed Einstein equations governing the d-dimensional 
on-vortex geometry therefore become 


— Kj Xp U — ^ t\ X v , 


(2.27) 
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of which maximal symmetry implies the only nontrivial combination is the trace 


R 


T^{d) 9 A R-IJ.U ~ tt/2 + 


d 


W 2 BW d 


BW'W 


d -1 


= -2TA, 


(2.28) 


and we use the explicit expression for R[d) i n terms of R and W. The components dictating the 
2-dimensional transverse geometry similarly are TZ mn = — K 2 X mn , which has two nontrivial 
components. One can be taken to be its trace 


ft(2) : = 9 


mn ft mn — R + d ( —— + 


IT" 

IT 


SW' 

BW 


= -K 2 X m m = -2 k 2 


6 ~ [1 -d 


X 


(2.29) 


and the other can be the difference between its two diagonal elements 

Q P P - g°0 = ft p p - n°e = -n 2 (t% - T°o) . 


(2.30) 


Writing out the curvature and stress energy shows this last equation becomes 


B (W'\ 
W \ld J 


2 

d 


k 2 Z. 


(2.31) 


Other useful combinations of Einstein equations 

Other linear combinations of the Einstein equations are not independent, but are sometimes 
more useful. The first is the {66) component of the trace-reversed Einstein equation lZ e g = 
—n 2 X e g which reads 


{B'W d y 

BW d 


Q-Z- 



= -2k 2 




(2.32) 


A second useful form is the (pp) Einstein equation, Q p p = —k 2 T p p , which is special in that 
all second derivatives with respect to p drop out. This leaves the following ‘constraint’ on the 
initial conditions for the integration in the p direction: 


B'W'\ R d{d- 1) fW'V 
BW ) + 2 W 2 + 2 [W) 


k 2 (z-x) 


K A gm L pot 



2.3 Scales and hierarchies 


(2.33) 


Before solving these field equations, we first briefly digress to summarize the relevant scales 
that appear in their solutions. The fundamental parameters of the problem are the gravita¬ 
tional constant, k; the gauge couplings, e (for Z M ) and g A (for A M )\ the scalar self-coupling, 
A, and the scalar vev v. These have the following engineering dimensions in powers of mass: 


[k] = 1-D/2, [ e ] = [g A \ = 2- D/2, [A] = 4 - D , and [v\ = D/2-1. (2.34) 
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To these must be added the dimensionless parameter, e, that measures the mixing strength 
for the two gauge fields. 

In terms of these we shall find that the energy density of the vortex is of order e 2 u 4 and 
this is localized within a region of order 


r v = 


ev 


(2.35) 


The effective energy-per-unit-area of the vortex is therefore of order e 2 v 4 r 2 = v 2 . These 
energies give rise to D-dimensional curvatures within the vortex of order 1 /L 2 = n 2 e 2 v 4: and 
integrated dimensional gravitational effects (like conical defect angles) of order k 2 v 2 . We work 
in a regime where kv <C 1 to ensure that the gravitational response to the energy density of 
the vortex is weak, and so, for example, defect angles are small and L v r v . 

By contrast, far from the vortex the curvature scale in the bulk turns out to be of order 
l/r| where 

A/V 

(2.36) 


r R = 


9a 


and J\f is a dimensionless measure of the total amount of A M flux that threads the compact 
transverse dimensions. Since our interest is in the regime where the vortex is much smaller 
than the transverse dimensions we throughout assume r v <C r B and so 


-^77 <C KV <C 1 
e,N 


(2.37) 


3. Isolated vortices 

We now describe some solutions to the above field equations, starting with the local properties 
of an isolated vortex within a much larger ambient bulk geometry. Our goal is to relate the 
properties of the vortex to the asymptotic behaviour of the bulk fields and their derivatives 
outside of (but near to) the vortex itself, with a view to using these as matching conditions 
when replacing the vortex with an effective codimension-2 localized object. These matching 
conditions are then used in later sections to see how a system of several vortices interact with 
one another within a compact transverse geometry. To this end we regard the field equations 
as to be integrated in the radial direction given a set of ‘initial conditions’ at the vortex 
centre. 

3.1 Vortex solutions 

For vortex solutions the vortex scalar vanishes at p = 0, and the vortex fields approach their 
vacuum values, if: —> v and 8 Z M —> 0, at large p. Because we work in the regime 1 these 

solutions closely resemble familiar Nielsen-Olesen solutions [3] in the absence of gravitational 
fields. Our analysis in this section reduces to that of [8] in the limit of no gauge mixing, 
e = 0, and a trivial bulk. 

8 In unitary gauge. 
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The asymptotic approach to the far-held vacuum values can be understood by linearizing 
the held equations about their vacuum conhgurations, writing ip = v + Sip and Zg = 0 + SZg. 
We find in this way that both Sip and SZ M describe massive particles, with respective masses 
given by 

m| = e 2 v 2 and m% = 2\v 2 . (3.1) 

From this we expect the approach to asymptopia to be exponentially fast over scales of 
order r z = m) 1 and rq, = m^ 1 . Indeed this expectation is borne out by explicit numerical 
evaluation. 

Notice the two vortex scales are identical, r v := r z = in the special BPS case, dehned 
by j3 = 1 where 

j3 := e 2 /2A , (3.2) 

and so the BPS case satishes e 2 = 2A. For convenience we also dehne (3 = e 2 /2A = (1 — e 2 )/3. 

Boundary conditions near the origin 

We start with a statement of the boundary conditions to be imposed at p = 0, which express 
that the transverse metric, g rnn , is locally hat and that all vectors (and so in particular the 
gradients of all scalars) must vanish there. For the metric functions we therefore impose the 
conditions 


W{ 0) = Wq , W'( 0) = 0 and B( 0) = 0 and B'( 0) = 1. (3.3) 

We can choose Wo = 1 by rescaling the d-dimensional coordinates, but this can only be 
done once so the change , AW, between the inside and the outside of the vortex (or between 
the centres of different vorticies) is a physical thing to be determined by the held equations. 
Similarly, for the vortex scalar we demand 

ip(0) = ip'(0) = 0, (3.4) 


or we could also trade one of these for the demand that ip —> v far from the vortex core. 
Nonsingularity of the bulk gauge field-strengths implies they must take the form 


A mn — f A Cr, 


Zmn — f z f-mn and SO A mn — /a C r 


(3.5) 


where e p g = = B is the volume form for the 2D metric g mn . Since A mn is nonsingular 

we know f A is regular at p = 0 and so because B(p) ~ p near p = 0 we see that A p g oc p near 
the origin. Consequently, in a gauge where A M dx M = A$(p) d 6 we should expect Ag = 0(p 2 ) 
near the origin. 

Naively, the same should be true for the vortex gauge helds A M and Z M , however the 
gauge transformation required to remove the phase everywhere from the order parameter 
T = ipe 1 ^ 1 (i.e. to reach unitary gauge) is singular at the origin, where T vanishes and so 
becomes ambiguous. Consequently in this gauge Zg (and so also A M ) does not vanish near 
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the origin like p 2 . Instead because in this gauge Z M —> 0 far from the vortex we see that flux 
quantization demands that 


27 m 


<f Z = 2vr f P ' 
J p=pv J o 


= § z (p < p v ) := f Z = 2n d pdpZg = 2n Zg(p v ) - Z e ( 0) = -2ttZo(0) , (3.6) 


where n is an integer, and we choose p = p v to be far enough from the vortex that 
there. We therefore ask Zg to satisfy the boundary condition: 

Tl TIE 

Zg( 0) = — and so therefore yWO) =-, 

e e 

where the second equality follows from Ag(0) = 0. 

Vortex solutions 

It is convenient to normalize the vortex fields 

Tl 

Zg = — P(p) and t/j = v F(p) 
e 


(3.7) 


(3.8) 


so that F = 1 corresponds to the vacuum value i/j = v, while the boundary conditions at 
p = 0 become 

F(0) = 0, P( 0) = 1; (3.9) 

the vacuum configuration in the far-field limit is 

F(oo) = 1, P( oo) = 0 . (3.10) 

In terms of P and F the Z M field equations boil down to 

1 


BW d 


W d P’\ e 2 v 2 F 2 P 


B 


B 2 


(3.11) 


while the ip equation reduces to 

— d ( BW d F'J = + An 2 F (F 2 - l) . (3.12) 

Although closed form solutions to these are not known, they are easily integrated numerically 
for given B and W, and the results agree with standard flat-space results when B = p and 
W = 1. See, for example, Fig. 1. 


BPS special case 

In the special case where W = 1 and e 2 = e 2 /(l — e 2 ) = 2A (and so /3 = 1), eqs. (3.11) and 
(3.12) are equivalent to the first-order equations, 9 


nP' 

Jb 


BF' = nFP and 


r 2 h 2 -i). 


(3.13) 












Figure 1: A comparison of BPS and non-BPS vortex profiles on a flat background for differing values 
of fj = e 2 /(2A). The (blue) profile vanishing at the origin is the scalar profile F and the (red) profile 
that decreases from the origin is the vector profile P. To find the profiles in flat space we set B = p and 
W = 1. The left plot uses $ = 1 and the right plot uses $ = 0.1, with this being the only parameter 
that controls vortex profiles in flat space. 


We show later that W = 1 also solves the Einstein equations when e 2 = 2A and so this choice 
provides a consistent solution to all the field equations in this case. 

When eqs. (3.13) and W = 1 hold, they also imply 

Ekin = - 2 (di ’) 2 = y i/j 2 Z m Z m = L gm , (3.14) 

and 

Lz := \ (1 - e 2 )Z mn Z mn = ^ 2 - v 2 ) 2 = V b , (3.15) 

which further imply that the vortex contributions to Z and X cancel out, 

£ = -hkin - -hgm = o and X loc = V b - L z = 0, (3.16) 

leaving only the bulk contribution to X: 


X = X B = A- L a . 


(3.17) 


As can be seen from eq. (2.31), it is the vanishing of Z that allows W = 1 to solve the Einstein 
equations. Finally, the vortex part of the action evaluates in this case to the simple result 


% : = 


v=§. 


d 2 yV^g 


Lqi + V b + L/2 


= 27r / dp B 


L<h + Vf) + Liz 


= nnv 


(3.18) 


9 The simplicity of these equations is understood in supersymmetric extensions of these models, since super- 
symmetry can require e 2 = 2A and the vortices in this case break only half of the theory’s supersymmetries. 
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P/r v 

Figure 2: A comparison of the profiles F and P for the vortex in flat space (dashed curves) and the 
full gravitating vortex solution (solid lines). For each case the (blue) profile that vanishes at the origin 
is the scalar profile F and the (red) profile that decreases from the origin is the vector profile P. The 
parameters used in the plot are d = 4, e = 0.3, (3 = 3, Q = 0.01 ev 2 , A = Q 2 / 2, kv = 0.6 and R = 0 
with the same values of /3 and e chosen for the non-gravitating solution. 


Bulk equations 


To obtain a full solution for a vortex coupled to gravity we must also solve the bulk field 
equations for W, B and A p g. The simplest of these to solve is the Maxwell equation, (2.16), 
whose solution is 



QB 
W d ’ 


(3.19) 


where Q is an integration constant. This enters into the Einstein equations, (2.28), (2.29) 
and (2.31), through the combination L A = /W d ) 2 . 

These can be numerically integrated out from p = 0, starting with the boundary condi¬ 
tions (3.3) (for which we choose Wq = 1), (3.4) and (3.6), provided that the curvature scalar, 
R, for the metric g pi , is also specified. Once this is done all held values and their derivatives 
are completely determined by the held equations for p > 0 and one such solution is shown in 
Fig. 3. As we shall see, many useful quantities far from the vortex depend only on certain 
integrals over the vortex prohles, rather than their detailed form. 
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Figure 3: These plots illustrate the bulk geometry for BPS vortices (/3 = 1) with parameters d = 4, 
e = 0, /3 = /3 = l , Q = 0.05 ev 2 , A = Q 2 /2 and kv = 0.3 (which also imply R = 0). In the top left 
plot, the solution for B is plotted (in blue) below the (red) metric function _B sp h e re of a sphere with 
radius r B = (200/3)r v . The presence of a vortex does not change the size of the bulk (since the full 
solution for B still vanishes at p = nr B ) and the metric function B is still approximately spherical with 
B ss 0.95 x £? S phere for these parameters. The top right plot shows that when /3 = 1 and A = Q 2 / 2, a 
constant warp factor solves the field equations. The bottom left plot shows that the derivative of the 
metric function B' ss 0.95 outside of the vortex core, at p > 4r v . The bottom right plot shows that 
B' ss —0.95 at the pole which lies opposite to the vortex core, indicating the presence of a conical 
singularity at that pole. 


3.2 Integral relations 

Our main interest in later sections is in how the vortices affect the bulk within which they 
reside, and this is governed by the boundary conditions they imply for the metric — i.e. on 
quantities like W, W ', B , B l and R — as well as for other bulk fields exterior to, but nearby, 
the vortex. In particular, simple integral expressions exist for derivatives of bulk fields — e.g. 
W' and B' — in this near-vortex limit, and we pause here to quote explicit expressions for 
these. 

For instance, imagine integrating the Einstein equation, (2.28), over the transverse di¬ 
mensions out to a proper distance p = p v — 0(r v ) outside of (but not too far from) the vortex 
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Figure 4: An illustration of the matching done at p — p v . The light grey surface is a cartoon of 
the bulk geometry. The bump on top of the surface represents the localized modifications to the 
approximately spherical bulk geometry that arise due to the vortex. The dark ring represents the 
circle at p = p v that lies sufficiently far outside the vortex that its fields are exponentially suppressed, 
but close enough to the vortex so that that its proper distance from the pole is still 0(r v ). 


(see Figure 4). This gives 


cLBW d d p In W 


P=Pv 


B ( W l 


p=p v 


= - — (2k z X + W- z R] 

2tt \ / v 


where we introduce the notation 


(0) v ■.= -^Jd 2 xV^gO = 2irJ* dp BW d O , 

x v 


(3.20) 


(3.21) 


and use the boundary condition W'( 0) = 0 at the vortex centre. This identifies explicitly the 
specific combination of vortex quantities relevant for specifying W' just outside the vortex. 

A second integral relation of this type starts instead with the (69) component of the 
trace-reversed Einstein equation: 1Z e q = — K 2 X e g , or (2.32), which integrates to give 

given the boundary condition B'W d —> 1 as p —> 0. This can be used to infer the implications 
of the vortex profiles on B' just outside the vortex. 

For many purposes our interest is in the order of magnitude of the integrals on the right- 
hand sides of expressions like (3.20) or (3.22) and these sometimes contain a surprise. In 
particular, naively one might think the integrals on the right-hand sides would generically 
be order v 2 and so would contribute at order k 2 v 2 to the quantities on the left-hand sides. 
Although this is true for g, the surprise is that the quantities (X) v and (Z) v can be much 



X — Z 


(3.22) 
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smaller than this, being suppressed by powers of r v /r B when the vortex is much smaller than 
the transverse space, r v <C r B , and this has important implications for how vortices influence 
their surroundings. 

One way to understand this suppression is to evaluate explicitly the suppressed quantities 
in the flat-space limit, where it can be shown (for instance) that the vortex solutions described 
above imply (T) v fl at = 0. Appendix A proves this as a general consequence of stress-energy 
conservation (or hydrostatic equilibrium) within the vortex, with the vortex dynamically 
adjusting to ensure it is true. (Alternatively, the vanishing of ( X) v on flat space can also be 
derived as a consequence of making the vortex action stationary with respect to rescalings 
of the size of the vortex.) More generally, for curved geometries we find numerically that in 
the generic situation when r v ~ r B all terms in (3.20) are similar in size and not particularly 
small, but this is no longer true once a hierarchy in scales exists between the size of the vortex 
and that of the transverse dimensions. In particular, as shown in Appendix A, for solutions 
where R is l/r| suppressed the vortex dynamically adjusts also to suppress ( X) v by powers 
of l/r B . 

The next sections provide several other ways to understand this suppression, associated 
with the constraints imposed by the Bianchi identities on the left-hand sides of near-vortex 
boundary conditions. 

3.3 Near-vortex asymptotics 

Because the vortex fields, 8^ = i/j — v and Z M , fall off exponentially they can be neglected 
to exponential accuracy ‘outside’ of the vortex; i.e. at distances p v r v ~ l/ev. The form 
for the metric functions B and W are then governed by the Einstein equations with only 
bulk-held stress-energy. This section describes the approximate form taken by these bulk 
solutions outside of the vortex sources, but not far outside (in units of the bulk curvature 
radius, say). 

We do so in two steps. We first solve for the bulk fields external to an isolated vortex 
in an infinite transverse space. We then find approximate asymptotic solutions for vortices 
sitting within compact spaces, under the assumption that the compact space is much larger 
than the transverse vortex size and that the region of interest for the solutions is very close 
to the vortex: r v <p«r B . 

Infinite transverse space 

We start with solutions where the space transverse to the vortex is not compact, since these 
should share many features of the bulk sufficiently close to a vortex residing within a large 
but finite transverse space. Concretely, the merit of seeking non-compact solutions is that 
the boundary conditions at infinity are fixed and determine many of the bulk integration 
constants. As seen in §4, compact spaces are more complicated from this point of view 
because these constants are instead set dynamically by the adjustment of the various vortices 
to each other’s presence. But those near-vortex boundary conditions that are dictated by 
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vortex properties should not care about distant details of the vortex environment, and so can 
be explored most simply within an isolated-vortex framework. 

To find isolated solutions we first write the Einstein equations in the exterior region 
p > p v where the vortex fields can be neglected: 


(W d B')' 

W d B 


fd-l 


Q 2 W~ 2d 



(3.23) 


and 

dB (jP) = 0 ’ ( 3 ‘ 24 ) 

and 

W~ 2 R + = k 2 [Q 2 W~ 2d - 2A) . (3.25) 

In this section (and only this section) we assume the transverse space does not close off, so 
B(p) >0 strictly for all values p > p v . 

Integrating (3.24) from p v to arbitrary p > p v gives 


W' W' 

—— = — = k 
B B v 


(3.26) 


where k is an integration constant and a v subscript indicates that the bulk field is evaluated 
at p = p v . Evaluating this at infinity tells us k = 0 if we demand W' vanishes there. More 
generally, if k ^ 0 and B monotonically increases then \ W\ must diverge at infinity, even if B 
(and so W') is bounded. Since B > 0 this excludes k < 0 since this would imply W vanishes 
at finite p > p v . If we also exclude W —> oo as p —> oo then we must have k = 0, for which 
integrating eq. (3.26) implies W = W v is constant everywhere outside the vortex. 

Using this result in eq. (3.23) gives 


B" 

13 


= —K 


d- 1 


q 2 w, 


-2d 


2A 
+ ~d 


=:Y d , 


(3.27) 


where the constancy of the right hand side (which we call Y d ) implies the transverse dimensions 
have constant curvature. Solving this for B in the region p > p v gives elementary solutions 
whose properties depend on the sign of : 

• Yd = —\/l 2 < 0: This implies B is a linear combination of sin (p/i) and cos (p/£) and 
so eventually passes through zero to pinch off at some r* > p v . This gives a compact 
transverse space, which we discard in this section. 

• Yd = +l/f 2 > 0 : This implies I? is a linear combination of sinh (p/£) and cosh (p/t) and 
so increases exponentially for large p. This corresponds to a vortex sitting within an 
infinite-volume transverse hyperbolic space with curvature radius l. 

• Yd = 0 : This forces B" = 0 which gives the flat solution B = B v + B[,(p — p v ). 
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A flat transverse space is found by tuning the bulk cosmological constant such that 
Yd = 0, and so 

A = — ^ (g? — 1) Q 2 W~ 2d < 0 . (3.28) 

Having A more negative than this gives a hyperbolic transverse space and more positive 
gives a compact transverse space. Evaluating (3.25) at the position p = p v , using (3.28) and 
constant W = W v then gives 

W~ 2 R = k 2 ( Q 2 W d - 2A) = dn 2 Q 2 W~ 2d = -2k 2 A ^ d ^ > 0, (3.29) 

which in our curvature conventions represents a strictly anti-de Sitter (AdS) geometry for the 
directions parallel to the vortex whenever the transverse directions are noncompact. 

As argued in more detail in §4.2, in general the 2D curvature scale, R = ±2/£ 2 , and the 
d D scale, R, are independent functions of the two dimensionful parameters: 1 /r\ oc k 2 A 
and 1/r 2 oc k 2 Q 2 . Of special interest is the one-parameter subspace of configurations for 
which either R or R vanish, and the above shows that the case R = 0 necessarily involves 
finite-volume transverse dimensions while flat transverse space {R = 0) implies an AdS on- 
vortex geometry, so the two subspaces intersect only as both r\ and r A tend to infinity (ie 
for A, Q —> 0). 

It is the constancy of W = W v in the bulk for isolated vortices that reflects something 
general about vortices: that W' —> 0 in the near-vortex limit. Indeed, although §4.2 gives 
explicit compact solutions with W' ^ 0 in the bulk, in all cases W' approaches zero in the 
immediate vicinity of the small source vortices. This carries implications for the integrated 
vortex stress-energy, such as ( X ) v . Using W' v = 0 in (3.20) allows us to write 

2k 2 {X) v = —R {W~ 2 ) v , (3.30) 

which is useful because it shows that {X) v is very generally suppressed by two powers of a 
curvature scale, being order p 2 /^ 2 <S 1 if R ~ l/£ 2 1 /p 2 . We expect this result also to 
hold for vortices situated within compact transverse dimensions. 

Asymptotic forms 

We next return to the case of real interest: small vortices situated within a much larger (but 
compact) transverse space. In general, the presence of a vortex introduces apparent singu¬ 
larities into the bulk geometry whose properties are dictated by those of the vortex. These 
singularities are only apparent because they are smoothed out once the interior structure of 
the vortex is included, since the geometry then responds to the stress-energy of the vortex 
interior. This section characterizes these singularities more precisely with a view to relating 
them to the properties of the source vortices. 

One way to characterize the position of the apparent singularity is to define it to occur 
at the point where the expression for Hb u ik(p) obtained using only the bulk field equations 
would vanish: i?buik(p*) = 0 (see Figure 5). Here p* is of order the vortex size, and need not 
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occur precisely at p = 0 (despite the boundary condition B(0) = 0 inside the vortex) because 
-Bbuik is found by solving only the bulk field equations without the vortex fields. 



Figure 5: A cartoon illustration of the definition of p*. The (blue) metric function B increases 
linearly away from the origin with unit slope B(p) « p. Outside of the vortex p ,> p v the solution is 
also linear in p but with B(p) ss ap. The straight (red) line extrapolates this exterior behaviour to 
the point, p = p+, where the external B would have vanished if the vortex had not intervened first. 


The nature of the singularity at p = p* is most simply described by expanding the bulk 
field equations in powers of proper distance, p = p — p*, away from the apparent singularity, 


W = W 0 
B = B 0 



(3.31) 


where r B is again a scale of order the bulk curvature scale. It is the leading powers, b and w, 
that describe potential singularity, and their form is constrained by the bulk field equations. 
In particular, as shown in Appendix C.l, the leading terms in the expansion of the Einstein 
equations around p = 0 imply that w and b satisfy the two Kasner conditions 10 [24]: 


dw + b = 1 and dw 2 + b 2 = 1. 


(3.32) 


The last of these in turn implies w and b must reside within the intervals 

|io| < —j=. and \b\ < 1. (3.33) 

Vd 

The Kasner solutions have precisely two solutions: either w = 0 and b = 1 (as is true for 
flat-space solutions) or w = 2/(1 + d) and b = (1 — d)/(l + d). Since we know that a non¬ 
gravitating vortex lives in a geometry with w = 0 and 6 = 1, this is also the root we must use 
in the weak-gravity limit ( kv ) 2 <S 1. This describes a conical singularity if B'(p = p*) ^ 1. 

1(, Our treatment here follows closely that of [25], which in turn is based on the classic BKL treatment of 
near-singularity time-dependence [26]. 
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Figure 6: Log-log plots of the the near vortex geometry for parameters d = 4, /? = 3, e = 0.3, 
Q = 0.01 ev 2 , A = Q 2 /2, kv = 0.6 and R = 0. The bulk in this case has a radius of r B = (500/3)r v . 
Outside of the vortex p > r v the geometry exhibits Kasner-like behaviour B' « a / 1 and W' ss 0. 


The field equations also dictate all but two of the remaining coefficients, B L and Wi, of 
the series solution. For instance eq. (2.31) applied outside the vortex implies W' = kB for 
constant k. This implies W\ = 0 and W 2 = \kar\ and so on, giving 



B = a p + ■ • • , 


P 2 + 


(3.34) 


where IT* = linip-^ W. For any such a singular point we therefore have the boundary 
conditions 

lim W' = 0 and lim B' =: a = const, (3.35) 

p->p* p->p* 

as is indeed found in detailed numerical integrations (see Figure 6). 

It is the slope B+ = a and LF* (where we affix W(0) = 1 within the vortex and so are not 
free to again choose IT* = 1 elsewhere) that convey the properties of the vortex to the bulk, 
and so should be governed by vortex properties, such as by boundary conditions like (3.20) 
or (3.22), rather than by bulk field equations. Notice that we expect both VF* — 1 and a — 1 
to be of order k 2 v 2 (see below) and so if W 2 = \kar 2 B is 0(1) then we expect k — 0(l/r|). 
This, in turn, implies that W' — 0(r v /r B ) at any near-vortex point of order r v away from p*. 
For r v r B we expect W to be approximately constant in the near-vortex region exterior 
to the vortex, up to 0{p 2 /r|) corrections. We also expect B' to be similarly constant, up to 
0(p/r B ) corrections. These expectations are verified by explicit numerical integrations of the 
vortex/bulk profiles, such as in Fig. 6. 

The explicit relation between a and vortex properties is set by near-vortex boundary 
conditions, such as (3.20) or (3.22). Using the series expansion to evaluate W and B at 
P = Pv, 

W = W v + Wy(p — p v ) -t- and B = B v + a (p - p v ) H-, (3.36) 
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where W v = IT* + ^ (p v — p*) 2 + •••, while W(, = kB v = ka(p v — p*) + ■ ■ ■ and so on. Inserting 
these into the left-hand side of eqs. (3.20) then gives 

dB v Wf~ 1 W[, = dkWt~ X c?pl + • • ■ = + W- 2 r) , (3.37) 

which confirms that the vortex adjusts to make the right-hand side CD(r 2 /r|). Similarly (3.22) 
becomes 

= + = (l~£) ’ (3 - 38) 

and so on. 


3.4 Effective description of a small vortex 

If the vortex is much smaller than the transverse space then most of the details of its structure 
should not be important when computing how it interacts with its environment. Its dynamics 
should be well described by an effective d-dimensional action that captures its transverse 
structure in a multipole expansion. 

The lowest-derivative ‘brane’ action of this type that depends on the nontrivial bulk fields 
outside the vortex is Sb = f d d x Cb with 


U = -V-l 


T _ — A , _1_ • • • 

± (£! C -^-pvAp I 


= -V=7 


p=Pb 


Til f mn A i 
-L i 2 fc ^mn W 


(3.39) 


J p=pb 


where 7 denotes the determinant of the induced metric on the d-dimensional world-volume of 
the vortex centre of mass (which in the coordinates used here is simply 7 ^ = g IJU evaluated at 
the brane position). The tensor A^ v \ p := ^ e p ,u\pmnA mn is proportional to the D-dimensional 
Hodge dual of the bulk field strength; a quantity that can be invariantly integrated over the 
d-dimensional world-volume of the codimension-2 vortex. All unwritten terms covered by the 
ellipses in (3.39) involve two or more derivatives. 

The dimensionful effective parameters T and ( respectively represent the vortex’s tension 
and localized flux, in a way we now make precise. To fix them in terms of the properties of 
the underlying vortex we perform a matching calculation; computing their effects on the bulk 
fields and comparing this to the parallel calculation using the full vortex solution. To do this 
we must be able to combine the d-dimensional action (3.39) with the D -dimensional action, 
S Bl for the bulk fields. 

To make this connection we promote (3.39) to a ZD-dimensional action by multiplying it 
by a ‘localization’ function, 5(y), writing the ZD-dimensional lagrangian density as 


£tot — kl B (g MN , A m ) + C b (g MN -, A M )S(y) ■ 


(3.40) 


Here C B is as given in (2.1) and 5(y) is a delta-function-like regularization function that 
has support only in a narrow region around the vortex position p = p b , normalized so that 
j v d 2 y 5(y) = 1. Although we can regard 5(y) as being independent of the d-dimensional 
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metric, g liu , and gauge field, A M , we cannot consider it to be independent of the transverse 
metric, g mn , because 6(y) must depend on the proper distance from the vortex. 

Much of the trick when matching with regularized delta-functions is to avoid questions 
that involve making assumptions about the detailed gWi-dependence of the brane action. This 
is most awkward when calculating the brane’s gravitational response, but we show below how 
to infer this response in a model-independent way that does not make ad-hoc assumptions 
about how 5{y) is regulated. 


Gauge-field matching 

We start with the determination of the coupling £ from the vortex’s gauge-field response. 

To determine £ we compute the contribution of S b to the gauge field equation, which 
becomes modified to 


dm 



6S b 
5 A n 


d„ 


V^9 A mn + (e r ' 


<%) 

y/9 2 


= o. 


(3.41) 


This has solution 


4 =9E 

P° yyd, 


C £p6» 


M 

y/92 


QB 

W d 


C<%), 


(3.42) 


where Q is an integration constant, and so — when integrated over a transverse volume, X v , 
completely containing the vortex — gives the flux 


<M^v) 


= I iA = Q I i2y {wi) 


C- 


(3.43) 


Comparing this to the vortex result in the full UV theory 

<MX V ) = <MC) -e* z (V) = Qj d 2 y ^ , (3.44) 

x v 


shows that ( is given at the classical level by 

2-rmE 

e 


(3.45) 


Notice that this argument does not make use of any detailed properties of S(y) beyond its 
normalization and independence of A m . 


Gauge-field back-reaction 

Before repeating this argument to match the tension, T, and determine the gravitational 
response, we first pause to draw attention to an important subtlety. The subtlety arises 
because the presence of localized flux causes the gauge field to back-react in a way that 
contributes to the localized energy density, in a manner similar to the way the classical 
Coulomb field back-reacts to renormalize the mass of a point charged particle. 
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To set up this discussion, notice that the effective lagrangian, (3.39), can be regarded as 
the macroscopic contribution of the vortex part of the lagrangian regarded as a function of 
applied fields A m and v . Consequently we expect the transverse average of (3.40) to give 
the same answer as the transverse average of the full lagrangian of the UV theory. Comparing 
the 7l m -dependent and -independent terms of this average then suggests the identifications 

T W d = (L kin + V b + L gm + L z ) v 

and ^ W d e mn A mn = = £ -(z™A m ^^ , (3.46) 

where W b = W(pb) is the warp factor evaluated at the brane position, and the factors W d 
come from the ratio of \/—7 /y/~9- 

Now comes the main point. The existence of the localized piece in the solution, (3.42), 
for A m has two related consequences in such a transverse average. 


• First, evaluating the localized-flux term at the solution to the A mn field equation, (3.42), 
shows that the localized component of A m renormalizes the tension, 


W b d 


T + | e mn A mn 


= W? 
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<%) 
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p=Pb 


(3.47) 


where this follows from taking 5(y) to be sufficiently peaked so that its integral can be 
treated like that of a Dirac delta-function. Notice that the last term in the last equality 
is singular as the vortex size goes to zero, requiring a regularization in order to be un¬ 
ambiguous. Such divergences are common for back-reacting objects with codimension-2 
or higher, and are ultimately dealt with by renormalizing the action (3.39) even at the 
classical level [27]. 

The ((-dependent part of this is to be compared with 

= S( Ls ) v , (3.48) 


which uses (2.18) and (3.19) to evaluate the integration over L m j x , and shows that the 
result agrees with (3.47), both on the value of the term linear in Q (once the matching 
value, (3.45), for £ is used) and by providing an explicit regularization of the singular 
0(e 2 ) term. 


The second way the localized term in (3.42) contributes is by introducing a localized 
contribution to the Maxwell action, L A , which was naively not part of the vortex 
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This exactly cancels the linear dependence on Q in (3.47), and partially cancels the 
localized renormalization of the tension. 


We see from this that the localized part of the gauge response to the brane action con¬ 
tributes a localized contribution to the bulk action (and energy density) that combines with 
the direct brane action in precisely the same way as happens microscopically from the mixing 
from A m to A m (see, for example, (2.19)). This suggests another useful notion of brane la- 
grangian, defined as the total localized contribution when Q is fixed (rather than A m ), leading 
to 


L b := TWg := (L kin + V b + L gm + L z ) = Wi 


c 2 

T - — 
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&(y) \ 

B K= Pb 


(3.50) 


We see that the tension renormalizations described above — associated with the [S(y)/B] Pb 
terms — are the macroscopic analogs of the renormalization e 2 —> e 2 = e 2 /(l — e 2 ) that occurs 
with the transition from L z to L z in the microscopic vortex picture. 

Whether L b or L b is of interest depends on the physical question being asked. L b arises 
in deriving the brane contribution to the A m field equations, as above. But because it is L b 
that contains all of the brane-localized contributions to the energy, it plays a more important 
role in the brane’s gravitational response (as we now explore in more detail). 


On-brane stress energy 

With the above definitions of L b and L b in hand we now turn to the determination of the 
brane’s local gravitational response. To determine the tension, T (or T), we compute the 
(/xzz) component of the Einstein equations (which we can do unambiguously because we know 
5(y) does not depend on g pv ). We can do so using either L b or L b to define the brane action. 
Using L b leads to the following stress energy 

r s >=£§ G£)=( r+ s £ra ”- 4 ™') % ^ ■ (3 - 5i) 

and so g becomes g = A + L A + g b with 

g b = < (t + ^ e mn A mn ^j ® . (3.52) 


Alternatively, using L b leads to the stress energy 
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(3.53) 


(3.54) 
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In either case the total energy density is the same, 


( Q ) v = {A + L A ) v + W b ^T+^e mn A. 


mn I — ( A + L a J + Wfr T , 
P=Pb 


(3.55) 


which is the analog of the microscopic statement (2.26) 

{o) v = ^A + L a + Akin + Lg m + Vb + L z + L m [ x ^ = ^A + L a + L\a n + T gm + V b + L z ^j . (3.56) 

The advantage of using (3.54) rather than (3.52) is that g b contains all of the brane-localized 
stress energy, unlike g b which misses the localized energy hidden in L A . 

IR metric boundary conditions 

A second important step in understanding the effective theory is to learn how the effective 
action modifies the held equations. So we restate here the general way of relating brane 
properties to near-brane derivatives of bulk fields [28]. The idea is to integrate the bulk held 
equations (including the brane sources) over a small region not much larger than (but totally 
including) the brane. For instance for a bulk scalar held, <h, coupled to a brane one might 
have the held equation 

□<F + J B + j b S{y) = 0 , (3.57) 

where J B is the contribution of bulk helds that remains smooth near the brane position and 
jb is the localized brane source. Integrating this over a tiny volume surrounding the brane 
and taking its size to zero — i.e. p v /r B —> 0 — then gives 

lim = 2vr lira B v Wf & = - lim ( J B + j b 6(y )\ = -j b {p = p b ) , (3.58) 

where the assumed smoothness of J B at the brane position ensures (J B ) V —> 0 in the limit 
p v 0. The equality of the second and last terms of this expression gives the desired relation 
between the near-brane derivative of <f> and the properties j b of the brane action. 

Applying this logic to the Einstein equations, integrating over a tiny volume, X v , com¬ 
pletely enclosing a vortex gives 
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(3.59) 


where we have split the action into an Einstein-Hilbert part S EH and a matter part S M . This 
matter part can be further divided into a piece that is smooth at the brane position 


S B = - J d D xsf^g (L a + A) , 

and one that contains all of the localized sources of stress energy. 
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As above, for a sufficiently small volume, X v , we need keep only the highest-derivative 
part of the Einstein-Hilbert term 11 , since the remainder vanishes on integration in the limit 
p v —> 0. The Sfj term also vanishes in this limit, by construction, so the result becomes 
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9 ik SS b \ 

$g jk / 


as p v —> 0 , 


(3.62) 


where i and j run over all coordinates except the radial direction, p, and K lJ is the extrinsic 
curvature tensor for the surfaces of constant p. To proceed, we assume that the derivative 
of the brane action is also localized such that its integral can be replaced with a quantity 
evaluated at the brane position 


/9^JiSb \ = f i 2 / g NP SS b \_f g NP 5S b \ 

\V^9 s 9 M p/ v J y V -J-g $g MP ) \V^Sgr) P=Pb 

X v 


(3.63) 


The b subscript in the functional derivative of the last term denotes that it is taken at the 
fixed point where 6(y) is localized, and so it contains no dependence on the bulk coordinates, 
and in particular no factors of 6(y). For example its pv components read 

(3.64) 

p=Pb 

However, at this point we remain agnostic about how to calculate the off-brane component 
5Sb/5gg0. Returning to the matching condition (3.62) we have the final result 
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which can be explicitly evaluated for the geometries of interest. 


(3.65) 


Brane stress-energies 

We now turn to the determination of the off-brane components of the brane stress-energy. 
We can learn these directly by computing the left hand side of (3.65) in the UV theory, before 
taking the limit p v —> 0. We will first do this very explicitly for the (pv) components of 
the brane stress-energy, and then proceed to deduce the off-brane components of the brane 
stress-energy. 

The (pv) stress-energy 

For the metric ansatz ds 2 = W 2 (p)g /iu daTMaT + dp 2 + B 2 (p)d0 2 , the extrinsic curvature 
evaluates to K t j = f g(-. This gives 

K> 1U = — g» v and K ee = B -g m . 

W B y 

11 Being careful to include the Gibbons-Hawking-York action [29] on the boundary. 


(3.66) 















Figure 7: A plot of defect angle matching in the region exterior but near to the vortex core. The 
solid (blue) lines represent the metric function W 4 B' and the dotted (red) lines represent 1 — k 2 T /2n 
computed independently for different values of £ = {—0.2, 0.2,0.4, 0.6} with the other parameters 
fixed at d = 4, f) = 3, Q = 1.25 x 10 _4 eu 2 , A = Q 2 / 2, kv = 0.5 and R = 0. This size of the defect 
angle B' v « a matches very well with 1 — n 2 T/2n at p = p v ~ 4r v . The solutions for W 4 B’ overlap 
perfectly when e = ±0.2, as indicated by the dashes in the line. This illustrates that the defect angle 
is controlled by T, and the linear dependence of the the defect angle on e is cancelled. 


The trace of the (pv) components of the condition (3.65) therefore evaluates to 
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for which the limit on the left-hand side can be evaluated using the limit B v —> 0 as p v —> 0. 
The result shows that it is the renormalized tension, T, that determines the defect angle just 
outside the vortex, 


k 2 W$T 
~2^ 


1 — a = 

This is the macroscopic analog of (3.22). 

The (98) stress-energy 

The (69) component of the metric matching condition, (3.65), evaluates to 
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but at first sight this is less useful because the unknown g mn dependence of 5(y) precludes 
evaluating its right-hand side. This problem can be side-stepped by using the constraint, 
eq. (2.33), evaluated at p = p v (just outside the brane or vortex) to evaluate W' v /W v = 
0(p v /r 2 B ) (and so also the left-hand side of (3.69)) in terms of the quantities B' v /B v = 1 / p v + 
• • •, R/Wy and X B . Once this is done we instead use the (99) matching condition to infer 
the (99) component of the vortex stress energy. 

Solving the constraint, (2.33), for W'/W at p v (just outside the vortex, where Z = 0 and 
X = X B = X B ) gives 
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(3.70) 


where the root is chosen such that W’ v /W v vanishes if both R and X B (p v ) vanish. With this 
expression we see that B v Wf ( W(,/W v ) —> 0 as p v 0, and so (3.69) then shows that 


(*"§) - 0 ' < 3 - 71 > 
x ' p=pb 

for any value of T (or T) and 

Notice that eq. (3.71) is precisely what is needed to ensure W' b —> 0 at the brane, as 
required by the Kasner equations (3.32) that govern the near-vortex limit of the bulk. Also 
notice that (3.71) would be counter-intuitive if instead one were to evaluate directly 5Sb/5ggg 
by assuming 5(y) was metric independent and using the explicit metrics that appear within 
e rnn A mn . What is missed by this type of naive calculation is the existence of the localized 
energy coming from the Maxwell action, L A , and its cancellation of the terms linear in £ when 
converting Sb to Sb- 

The (pp) stress-energy 

Although the (pp) component of the extrinsic curvature tensor is not strictly well-defined, we 
can still consider the (pp) components of the boundary condition in the following form 
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By definition, we have 
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As noted in (2.33), this component of the Einstein tensor is contains only first derivatives of 
the metric field. It follows that 
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since metric functions and their first derivatives are assumed to be smooth. In this simple 
way, we once again use the Hamiltonian constraint to conclude that the off-brane component 
of the brane stress energy is vanishing 


( g pp ss b 

\V=§stf 


= o. 

p=pb 


(3.75) 


So both off-brane components of brane stress-energy vanish in the limit p v —> 0, and from 
this we also infer that their sums and differences also vanish: 


X b = Z b = 0 . (3.76) 

These results are the analog for the effective theory of the KK-suppression of (X) v in the UV 
theory once r v <r B . As a consequence in the effective theory 

Z = 0 and X = X B = A- L A . (3.77) 

4. Compactification and interacting vortices 

We next turn to how several small vortices interact with one another and with their environ¬ 
ment. In particular, if the flux in the transverse dimensions does not fall off quickly enough its 
gravitational held eventually dominates and drives B(p) to zero for positive p , thereby pinch¬ 
ing off and compactifying the two transverse dimensions. We explore in detail the situation 
of two small vortices situated at opposite sides of such a compact space. 

For this part of the discussion it is more convenient to use the effective description of 
the vortices as codimension-2 branes than to delve into their detailed vortex substructure, 
though we do this as well to see how the effective description captures the full theory’s low- 
energy behaviour. As we saw above, in the effective limit the vortex properties are encoded 
in the near-brane derivatives of the bulk fields (through the defect angle and localized flux). 
So to discuss brane interactions it is useful to start with the general solution to the bulk 
held equations outside the vortices, since it is the trading of the integration constants of 
this solution for the near-brane boundary conditions that expresses how brane properties 
back-react onto their environs. 

4.1 Integral relations 

Before delving into explicit solutions to the bulk held equations, it is worth hrst recording 
some exact results that can be obtained by applying the integrated forms of the held equations 
to the entire transverse space, and not just to a small region encompassing each vortex. In 
the UV theory these integrals simplify because all helds are everywhere smooth and so the 
integral over a total derivative vanishes. The same need not be true for the effective theory 
with point brane sources, since in principle helds can diverge at the brane locations. However 
we can then ask how the UV integral relations arise in the effective theory. 
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For instance if eq. (3.20) is integrated over the entire compact transverse space then its 
left-hand side integrates to zero, leaving the following exact relation for R : 

0 = (2 k 2 x) +(w~ 2 ) R = 2k 2 fWtot + AV (4.1) 

\ / tot \ / tot y 

Here the last equality uses the relation between k 2 and its d-dimensional counterpart. This 
shows that it is (A’)tot that ultimately determines the value of the on-brane curvature. 

Eq. (4.1) is particularly powerful in the effective theory, for which we have seen that the 

branes satisfy Xb = 0 and so X = X B = A — L A . In this case (4.1) shows us that it is really 

only through 

{La )tot = 27 tQ J dp (4-2) 

that the brane properties determine the on-brane curvature, as they modify the functional 
form of B and W d through boundary conditions, and Q through flux quantization. 

A second exact integral relation comes from integrating the {99) component of the trace- 
reversed Einstein equation, eq. (2.32), over the transverse dimensions. Again the left-hand 
side integrates to zero leaving the constraint 

(4 - 3 > 

Combining this with (4.1) then implies 

R = ~ (d^) K 2 d ( g ~ Z )tot' (44) 

Notice that in d dimensions Einstein’s equations with a cosmological constant, V e R, have 
the form 

Rjiv — T) = K d ^eff 1 (4-5) 

and so the scalar curvature satisfies 

k = - (5^2) ■ (4 6 > 

Comparing this with eq. (4.4) then gives a general expression for the effective d-dimensional 
cosmological constant 

V eS = ( e -z) . (4.7) 

\ / tot 

4.2 General static bulk solutions 

This section presents the general bulk solutions for two branes. We start with the simple 
rugby-ball geometries that interpolate between two branes sourcing identical defect angles 
and then continue to the general case of two different branes. The solutions we find are 
all static - actually maximally symmetric in the d Lorentzian on-brane directions - and 
symmetric under axial rotations about the two brane sources. 
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Rather than rewriting all of the field equations for the bulk region away from the branes, 
we note that these are easily obtained from the field equations of previous sections in the 
special case that Z M = 0 and ip = v. Notice that Z M = 0 and ip = v already solve the Z M 
and ip field equations, so it is only the others that need solutions, which must be the case 
since we have replaced the vortex degrees of freedom with an effective brane description. 

These choices imply 

1 / Q \ 2 

hkin — Lty — Vb — hgm — L z — Tmix — 0 and so Tgg e — L A — L A — 2 ( yyd J ' (4-8) 

As a consequence of these we know 

Z = 0 , X = X B = A — L a and g = g B = A + L A . (4-9) 


Rugby-ball geometries 

Because Z = 0 the solutions to the field equations can be (but need not be) locally maximally 
symmetric in the transverse 2 dimensions, rather than just axially symmetric there. For such 
solutions W' must vanish and so the geometry is completely described by the constant scalar 
curvatures, R and R. The transverse dimensions are locally spherical, but with a defect angle 
at both poles corresponding to the removal of a wedge of the sphere. 

Explicitly, we have B = a£sin(p/£), and the polar defect angle has size 6 = 27r(l — a). 
The sphere’s curvature and volume are 


2£// 2 /*7r£ 

R = —— = — —j- and tb :=2tt dp B = 47tqT 2 , 

B t Jo 


(4.10) 


where £ is the ‘radius’ of the sphere. The relevant bulk field equations are the two Einstein 
equations 
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with Q fixed by flux quantization to be 

Q _ M 
W d ~ 2g A a£ 2 


where N := N — n to t£ ^ 


(4.11) 

(4.12) 

(4.13) 


where ntot = + n- is the sum of the flux quantum for each vortex. 12 

As shown in Appendix C.2, the stable solution to these equations has compact transverse 
dimensions with radius 



(4.14) 


12 We take for simplicity the gauge coupling of the two vortices to be equal. See Appendix B for a discussion 
of flux quantization for the Z M and A M fields. 
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where the two intrinsic length-scales of the problem are defined by 


r A := 


4a 2 A 


and 



(4.15) 


Clearly £ ~ r A /yJ2 when r\ S> r A and increases to £ = r A when r A = r A . It is here that 
we first see why it is the combination Mn/g A that sets the size of the extra dimensions. No 
solutions of the type we seek exist at all unless r A > r A , which requires 
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Finally, the on-brane curvature is 
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(4.16) 


(4.17) 

(4.18) 

(4.19) 


The on-brane curvature passes through zero when A is adjusted to satisfy r 2 A /r\ = 4 (d — 1) /d 2 
(which is < 1 for d > 2), and £ 2 = := (2 /d)r\. 


Geometries for general brane pairs 

Explicit closed-form solutions are also known where the branes at either end of the space have 
different properties. The difference between the two branes induces nontrivial warping and 
thereby breaks the maximal 2D symmetry of the transverse dimensions down to simple axial 
rotations. 

The resulting solutions can be found by double Wick-rotating a charged black hole solu¬ 
tion in D dimensions [30, 31], leading to the metric 


( di9 2 


dfig = W 2 (d) g^ dx^dx" + r% + a o K ($) sin 2 ddd 2 

= W 2 (d) g^y dx^dx” + r 2 (d) (dd 2 + a 2 (d) sin 2 d d 0 2 ^j , 


where 


(4.20) 


W{d) := Wg^l + g cosi?^ , r(d) := ~^ = = and a(d) := ao K(d ), (4-21) 
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where g is an integration constant and r^ 2 := 2/t 2 A = (d/2)r^ 2 . Notice that r 2 (i?)a($) = 
7 q ao for all $, and the vanishing of ggg implies the ‘radial’ coordinate lies within the range 
i9 + := 0 < i? < := 7 r. The geometry at the endpoints has a defect angle given by 

a± = a(i?±) and the derivative of the warp factor vanishes at both ends: dW/dJ) —> 0 as 
i? —> d± (as required by the general Kasner arguments of earlier sections). In these coordinates 
the Maxwell field solves \J—g A® e = Q , which implies 


A$e 


Q Tq ao sin $ 

wRj) 


(4.22) 


Other properties of this metric — including the explicit form for the function K('&) — are 
given in Appendix C.2. 

All told, the solution is characterized by three independent integration constants, in terms 
of which all other quantities (such as R ) can be computed. These constant can be taken to 
be Q as well as an independent defect angle, a + and a_, at each of the two poles. These 
three constants are themselves determined in terms of the source brane properties through 
the near-brane boundary conditions and the flux-quantization condition 
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(4.23) 


where, as before, M 


N—ntot£g A /e represents the vortex-modified flux-quantization integer. 


Near rugby-ball limit 

Although the general expressions are more cumbersome, it is possible to give simple formulae 
for physical quantities in terms of Q and a± in the special case where the geometry is not 
too different than a rugby ball. Because nonzero g quantifies the deviation from a rugby-ball 
solution, in this regime we may expand in powers of g. In this section we quote explicit 
expressions that hold at low order in this expansion. 

In the rugby-ball limit the functions r(i9), a(i?) and W(tJ) degenerate to constants, with 
W{d) = Wo and r(i?) = l given explicitly in terms of Tq = (2 /d)r\ and R by eq. (4.17). Since 
r(i?) 2 a(i9) = Tq ao this implies ao is related to the limiting rugby-ball defect angle, a, by 


a = ao 


1 + 


d- 1 



(4.24) 


It happens that to linear order in g the geometry near each pole takes the form 


ds^_ ~ Wq (1 ± 2 g)ds\ + 


dr? 2 + <4(# - *?±) 2 d6 2 


(4.25) 


where 


AW -=W + -W^~2Wog + 0{g 2 ) 

4 e (l ± C H g ) + O [ 77(0 - tf±) 2 , g 2 ] , (4.26) 

a± ~ a(l^C H g^j + O [g(J} -‘d±) 2 ,g 2 ] . (4.27) 
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with 


C H : = 


d — | + (d — 1 )Rvq 


(4.28) 


1 — (d — l)Rr'^/d 

This shows that the apparent rugby-ball radius and defect angle seen by a near-brane observer 
at each pole begins to differ for each brane at linear order in g. 

To use these expressions to determine quantities in terms only of Q and a± requires 
knowledge of R, and the field equations imply this is given for small g by 


R =-FT 
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To complete the story we solve for in terms of a± using 

a_ — a_|_ 
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— 2 C H g , 


with a — ^(aq_ + a_), and use this to evaluate all other quantities. 

For small g the flux-quantization condition also simplifies, becoming 
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4.3 Relating bulk to vortex properties 
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1 + d(d+v 2+ < 
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(4.29) 


(4.30) 


(4.31) 


We see that the bulk solutions are determined by three parameters, ot± and Q. Earlier sections 
also show how these are related to the physical properties of the two source branes, since the 
near-brane defects are related to the renormalized brane tensions by 


1 — a± 


k 2 W%T± 
2tt 


(4.32) 


and Q is determined in terms of brane properties by flux quantization, (4.23) (or, for small 
V, (4-31)). 

Parameter counting 


An important question is to count parameters to see if there are enough integration constants 
in the bulk solutions to interpolate between arbitrary choices for the two vortex sources. 

In total the source branes are characterized by a total of four physical choices: their 
tensions (i.e. defect angles) and localized flux quanta, n±, to which we must add the overall 
flux quantum choice, N, for the bulk. But varying these only sweeps out a three-parameter 
set of bulk configurations because the flux choices ( n± and N ) only appear in the bulk 
geometry through flux quantization, and so only through the one combination, Af = N — 
e(n + + n-)(g A /e), that fixes Q. (Although they do not affect the geometry independent of 
Af, the n± do govern the Bohm-Aharonov phase acquired by test particles that move about 
the source vortices.) 

Consequently the three free constants — Q and a± — are sufficient to describe the static 
gravitational held set up by any pair of vortices, and once the brane properties (and N ) 
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are specified then all geometrical properties are completely fixed. The rugby ball geometries 
correspond to the special cases where T + = T_. 

This point of view, where the bulk dynamically relaxes in response to the presence of 
two brane sources, is complimentary to our earlier perspective which regarded integrating the 
field equations as an ‘evolution’ in the radial direction (and so for which initial conditions at 
one brane completely specify the whole geometry — and by so doing also fix the properties 
of the antipodal brane). They are equivalent because in the evolutionary point of view two of 
the integration constants to be chosen were Q and R, which are completely arbitrary from the 
perspective of any one brane. Their choices dictate the form of the interpolating geometry 
and correspond to the two-parameter family of branes (labeled by A f and a ) that could have 
been chosen to sit at the antipodal position. 

On-brane curvature response 

Of particular interest is how the on-brane curvature, R, responds to different choices for brane 
properties. In general this is given by (4.4), in which we use the brane-vortex matching results 
- (3.55) and (3.76) — appropriate when the vortex size is negligibly small compared with 
the transverse KK scale, Z b = X b = 0 and g b = W b T b 8(y)/ ensuring that 

(^)tot — 0 ; ('T)tot — (A — L a ) t 0 t and (g)tot = W :b T b + (A + L A ) to t ■ (4.33) 

b 

With these results (4.4) shows R takes the value that would be expected in d dimensions in 
the presence of a cosmological constant of size V e g = (g)tot = (1 — 2 /d) (X) tot , and so 

Kff = ^ Wf T b +(A + La) tot = (l - ^ (A — L a ) to t, (4.34) 

In general R is not small. Since all quantities in g are positive (except perhaps for A), the 
resulting geometry is de Sitter-like unless cancelled by sufficiently negative A. Notice also 
that the second equality implies 

J2 W b d T b = -- d (A) tot — 2^1 — £) <L A ) tot , (4.35) 

is always true. This states that for codimension-2 systems like this the ‘probe’ approximation 
is never a good one: that is, it is never a good approximation to neglect the bulk response 
(the right-hand side) relative to the source tensions (the left-hand side) themselves. 

Near-flat response 

Of particular interest are near-flat solutions for which A is initially adjusted to cancel the 
rest of (g)tot, after which brane properties are varied (without again readjusting A). One 
can ask how R responds to this variation. To determine this response explicitly we use the 
near-rugby solution considered above, in the case where the unperturbed flat geometry is a 
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rugby ball and for which the brane parameters are independently tweaked. To this end we 
take the initial unperturbed configuration to satisfy Wq = 1 and 


O 2 

A = Aq := and rjo = 0 


Rq — 0 


(4.36) 


and then introduce small perturbations through 5a and 5J\f. 
diately that 


5R 


l 5 ® 

H) Qo 


+ o( v 2 ) 


From eq. (4.29), we see imme- 


(4.37) 


and — from the flux quantization condition in eq. (4.31) - 
bations are 


5Q 

Qo 


5M_ 

■A/q 


5a® 

ao 


+ 0(rj 2 ). 


we see that the leading pertur- 


(4.38) 


Lastly, since it is a = ^(a+ + a_) (and not «o) that is determined by T±, we must use the 
relation — eq. (4.24) — to write 


5a o 
a o 




5R + 0(rf). 


Combining these formulae gives 


d- 2 





to leading order, and so given perturbations of the form 


5a± = 


k 2 5T± 

2ir 


5U 


fce± 


(4.39) 


(4.40) 


(4.41) 


the corresponding change in R has the form 


d- 2 
d 




o 


b=± 


K 2 5T b 

2na 




(4.42) 


Comparing with (4.6) — and using k^ 2 = Aitar^/ k 2 and k 2 QqTq = 1 for the unperturbed flat 
rugby-ball geometry — then shows that this curvature is what would have arisen from the 
d-dimensional vacuum energy 

= ( stb ~ • ( 4 - 43 ) 

We see from this that when 5J\f = 0 the curvature obtained is precisely what would be 
expected in d dimensions given the energy change 
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5. Discussion 


In this paper we investigated the gravitational properties of branes that carry localized flux 
of a bulk held, or BLF branes. As noted in the introduction, the treatment of a gravitating 
BLF branes is not straightforward because the delta-like function used to represent their 
localization must depend on the proper distance away from the brane. Because of their 
particularly simple structure, this is not a problem for branes described only by their tension 
oc T. However, the presence of metric factors in the BLF term oc e rnn A mn complicates any 
calculation of transverse components of the brane’s stress energy. 

We resolved this ambiguity by constructing an explicit UV completion of BLF branes 
using Nielsen-Oleson vortices whose gauge sector mixes kinetically with a bulk gauge held. 
The gauge kinetic mixing, which is controlled by a dimensionless parameter e, endows the 
bulk held with a non-zero hux in the localized region, even in the limit that this region is taken 
to be vanishingly small. This allows the UV theory to capture the effects of brane-localized 
hux. 

The main result is that, in the UV picture, the gauge kinetic mixing can be diagonalized 
resulting in variables that clearly separate the localized sources from the bulk sources. In 
the diagonal basis, the energy associated with localized hux is always cancelled, and the 
canonical vortex gauge coupling is renormalized: e 2 = e 2 /(l — e 2 ). This allows us to identify 
the renormalized vortex tension as the quantity that controls the size of the defect angle in 
the geometry exterior to the vortex. We also hnd that the vortex relaxes to ensure that the 
average of the localized contributions to the transverse stress energy are controlled by the 
ratio between the size of the vortex and the characteristic bulk length scale r v /r B . 

This informs our treatment of the IR theory with branes. We hnd that the delta-function 
treatment of the brane is particularly useful for calculating the hux of the bulk held, includ¬ 
ing its localized contributions, and a delta-function shift in the bulk gauge held strength can 
diagonalize the brane-localized hux term. This change of variables endows the action with a 
divergent term that we can interpret as a renormalization of the brane tension, in analogy 
with the e —> e renormalization of the gauge coupling. We also show that the transverse com¬ 
ponents of the brane stress energy must vanish without explicitly calculating them. Rather, 
we use the Hamiltonian constraint and energy conservation to relate these stress energies to 
quantities which vanish as r v /r B —> 0, thereby circumventing any ambiguity in the metric 
dependence of the corresponding brane interactions. 

The techniques we employ here should be relevant to other brane couplings that contain 
metric factors. For example, there is a codimension-A: analogue of the BLF term that involves 
the Hodge dual of an fc-form. Of particular interest is the case k = 1 where the brane can 
couple to the derivative of a bulk scalar field (f> as follows Sb oc f *d</>, or a bulk gauge field A 
as follows Sb oc f -kA. We have also provided an explicit regularization of a <5(0) divergence. 
These are commonplace in treatments of brane physics, and usually deemed problematic. 
However, there is likely a similar renormalization story in these other cases. 

Lastly, we plan [21] to also apply these techniques to a supersymmetric brane-world 
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models that aim to tackle the cosmological constant problem [20]. The back-reaction of branes 
is a crucial ingredient of such models, and understanding the system in greater detail with 
an explicit UV completion will put these models on firmer ground and hopefully shed light 
on new angles from which to attack the CC problem. 
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A. Stress-energy conservation 


The matter held equations always guarantee the matter stress energy is covariantly conserved, 
V m T mjv = 0. For the geometries of interest this has one nontrivial component, \7 M T Mp = 0, 
which implies 

[BW d TP„y = BW d T°o + ^ T>\\ . (A.l) 

A useful way to rewrite this multiplies by B and adds BB'W d T p p to both sides, so 


B 2 W d T p p\ = BW d 


B' ( T\ + T p p ) + ^ T\ 


(A.2) 


or 


/ 


V^gB(z-x) 


-V=g 


2 B'X + 


dBW' 

W 



When applied to a vortex on hat space — for which W = B' = 1 and the constraint 
implies Z — X = 0 outside the vortex — integrating eq. (A.3) over the vortex reduces 
simple statement 


(A.3) 

(2.33) 
to the 


Wvlflat “ ( Lpot L SS e )vlflat “ °’ 


(A.4) 
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a result that may also be derived as the vortex equation of motion corresponding to extrem- 
izing the flat-space action against rigid rescalings. 


B. Flux quantization 


For compact transverse dimensions the underlying transverse geometry of interest has the 
topology of a sphere, leading to two flux-quantization conditions; one for each of the U{ 1) 
gauge fields. 

Our interest is in vortices that are much smaller than the size of the transverse space. 
In this case we take the complex vortex field and the gauge field Z M at the equator to be 
gauge-equivalent to their vacuum values, \F = ve ieu and Z M = d M u. Here single-valuedness 
of \F± on both patches, V±, implies 

2irn± 


u±(6 + 2i r) = u±(9) + 


(B.l) 


for some integers n±. The choices for n± can be chosen differently because they differ by a 
gauge transformation, g = e l6UJ , whose single-valuedness implies lo(9 + 2ir) = ui(9) + 2itN z , 
provided N z = n + + n_ =: n to t- 

The total Z-flux through each hemisphere is related to the integral of Zq around the 
equator by 

2 '"^. (B.2) 


(B.3) 


<f> z ± = dZ = ® Z = 

Jv± Jdv± 


and so the total Z-flux through the sphere is 


<£„ = 


2irn 


tot 


For the A M gauge field we imagine test charges situated far from the vortices that couple 
to A m and carry charge g A . The action for such a charge probe is 


Sprobe — Qa / A , 


(B.4) 


exp 

i 

1 ’ ^ 

= exp 

ig A / d A 

= exp 

i- 

•<s> 

i 


L Je q J 


[ Jv+ \ 


L JV- J 


where the integration is along the world-line of the charge. Moving such a charge around the 
equator far from the vortex contributes an amount 


(B.5) 


to the path integral, where the two equalities rewrite the integral using Stoke’s theorem and 
the observation that the equator can be considered to be the boundary of either hemisphere 
(with the sign keeping track of the orientation of the boundary). In order for this phase to 
be single-valued in the path integral we must therefore demand the fluxes, = fp dH, 
satisfy 

$ A := $ A+ + $. 4 - = , (B.6) 


9a 


for some integer N. 
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Bulk vs localized flux 

Suppose now we take a test charge that starts life coupled only to A M and move it around the 
vortex, keeping always far enough from the vortex that the Z M magnetic field is negligible. 
Then we define the flux seen by this charge by 


$. 4 pf v ) := (p A= I dA= I (dA-edZj 

27TO 


ax, x, x v 

p " f QB ^ ne 


= 2vr 


Uo 


dp 


W d J 


(B.7) 


where the second-last equality uses flux quantization (for integer to) for the vortex solution 
for Z localized well within region V : 


$z(X v ) := / d Z = - 


2irn 


(B. 8 ) 


x v 


and the sign on the far right-hand side is chosen for later convenience. The approximate 
equality in (B.7) drops the order ( p v /t ) 2 contribution of the A flux over the vortex volume 
relative to the localized Z flux. 

What is important about (B.7) is that the gauge-field mixing implies that the test charge 
now sees a vortex-localized component due to the appearance of the Z term. It is in this 
sense that our system provides a UV completion for branes carrying nonzero brane-localized 
flux. 

On general grounds the flux of A is also quantized, and this fixes the value of Q. That 
is, if the integration is performed over the entire transverse dimensions it must satisfy 

, . 27rN , 

$4 (tot) = -, (B.9) 

9a 


where N is an integer and g A is the gauge coupling of the field A M to its test charge. Con¬ 
sequently, the presence of the brane-localized flux modifies what flux quantization demands 
for Q: 



(B.10) 


where p = 1* denotes the proper distance between the branes (defined by the two places 
where B vanishes: B( 0) = -E>(7*) = 0), and ntot = n + + n - is th e sum of the flux quanta for 
the vortices localized at each of these positions. 


C. Solutions 

This appendix describes more details about the approximate and exact solutions described 
in the main text. 
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C.l Approximate near-vortex solutions 

For the purposes of matching the bulk integration constants to the vortex properties we are 
most interested in the form of the solutions very near to, but outside of, the vortex sources. 
We start by recapping the form of the bulk solutions very close, but outside of, a small vortex. 


Asymptotic forms 

Near the branes it is possible to expand the solutions in powers of p/r B , where p denotes 
proper distance in the bulk geometry from the vortex. Writing, as before, the metric in the 
form 

ds 2 = W 2 {p) dx^dx" + dp 2 + B 2 (p) dd 2 , (C.l) 

we seek near-vortex solutions to the Einstein equations of the form 


W = W 0 
B = B 0 



(C.2) 


and so on. The special case of flat space in polar coordinates corresponds to w = 0 and 6 = 1, 
without the need for higher powers of p/r B . 

The leading powers, w and b, are constrained by the leading terms in the expansion of 
the field equations around the vortex position, p = 0. The source terms on the RHS of the 
Einstein equations in the bulk involve A and L A = \{Q/W d ) 2 , which vary respectively like p° 
and ( p/r B )~ 2dw as p —>• 0. By comparison, as p 0 the curvature on the LHS of the Einstein 
equation are 


R 


n ( d ) W 2 


—r = d 


,, lW W'\ 2 f W" B'W' 
^ ~~ ^ ~W ) + (“W + BW 


= d{(d-l)[- p ) + 


iv 


= dw 


dw + b — 1 
P 2 


w{w — 1) bw 

J 2 + 


P 


1 + 0 ( ^ 

r B 


1 + 0 ( dh 

r B 


(C.3) 


Assuming w < 1 — so that R/W 2 oc ( p/r B )~ 2w is subdominant to the 1/p 2 term explicitly 
displayed (a result justified below) — we see that the (pv) Einstein equation implies w(dw + 
b — 1) = 0. Similarly, 


D// 

n e e = ^ + d 

implies b(dw + b — 1) = 0, and 


B'W'\ _ b(dw + b-l) 
BW 


1 + o(^ 

r B 


vp _B^_.JW^\_Kb-l) + dw{w-\) 
K p ~ B + d l W 1 “ 


1 + 0 ( ^ 

r B 


(C.4) 


(C.5) 
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Besides the trivial special case {w = 6 = 0) we see that the vanishing of the 1/p 2 terms 
in the held equations implies the following two Kasner conditions: 


dw + 6 = 1, 


(C.6) 

and 



dw 2 + 6 2 = 1 


(C.7) 

The last of these in turn implies w and 6 must reside within the intervals 


ke < and 

Vd 


(C.8) 

which shows in particular why l/W 2 oc ( p/r B )~ 2w 

is less singular than 1/p 2 , 

as assumed 


above. The Kasner conditions, eqs. (C.6) and (C.7) have precisely two solutions: either 
w = 0 and 6=1 (as is true for the rugby-ball solutions described above) or dw = 1 and 6 = 0. 

C.2 Exact Solutions 

This section explores some properties of the solutions to the vortex-bulk held equations. 


Rugby-ball geometries 

We next describe some details associated with the rugby-ball geometries, for which solutions 
to the held equations with = v and Z M = 0 are sought of the form B = a£sm(p/£), with 
W constant. The transverse curvature and volume are 


R = 


2 B' 1 

If 


2 

'F 


and 


PIT0. 

V 2 :=2ir dp B = 4ira£ 2 , 

Jo 


(C.9) 


where l is the ‘radius’ of the sphere. The relevant bulk held equations are the two Einstein 
equations 

.2 


and 


R = -2k z 


R ^ 2k 


A — — ( -Q— 
2 \ W d 


2A 

T + 11 - 


d \w d 


with Q hxed by hux quantization to be 

Q _ M 
W d 2 g A a£ 2 


where AT := N — ntot£ (J^j 


(C.10) 


(C.ll) 


(C.12) 


Using (C.12) in (4.12) allows £ 2 to be solved as a function of brane properties (which 
enter through the defect-angle parameter a). Dehning 


o d 

r\ : = 


4ac 2 A 


and 


t a(oi) ■■= - 1 - -7 


1 


Mk 


dj \g A a 


(C.13) 
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we find 

_ y? + ~T = 0 or > equivalently £ 4 - 2r\i 2 + r\r\ = 0 , (C.14) 

Hr r\ 

which has solutions 

= ^1 ± jj 1 - ^ or ’ equivalently 4 = r| ^1 T ^1 - ^ ■ (C.15) 

Clearly tL ~ \/2 ta and t?+ ~ r A /\J~2 when ?’a 2> r A . As ta decreases from infinity also 
decreases and £+ increases until they meet at t+ = £_ = r A = ?’a when ta = r A . No solutions 
of the type we seek exist if ?’a < r A , and so the existence of solutions requires we choose 


d-Wafo>,2 

“ 2 W/cV 

Finally, the corresponding on-brane curvature becomes 


(C.16) 


R± = 


d 


+ 


d 


2 r\ 2{d-l) V4 

which shows 


+ 


2(d- 1)4 Vrf- 1 


-w l± 4/1 — -4 


i?_ ~ — 


d(2d —1)1 1 

L2(rf-1) J 4’ 


and Ab ~ 


2d 

d — l ) r A 


2 


if r A > r A . 


(C.17) 


(C.18) 


These two roots converge towards one another as ?’a r A , eventually converging at 


R = - 


d-2\ 2 

~dZT\ ) ^2 when r A = r A ■ 


(C.19) 


Although AL remains negative (de Sitter space) between these two limits, R + changes sign. 
It passes through zero when 


d 2 


(which is < 1 for d > 2). 


(C.20) 


The origin of the two roots can be understood by thinking of the lagrangian as a function 
of l in order to understand the potential that i is minimizing. Since L tot = L EH + L A + A, 
and since 77( 2 ) = R = — 2/i 2 and L A oc r 2 A /i 4 , and since oc l 2 , we can see that in the 4D 
Einstein frame we have a potential for £ that involves three terms: 


V eS (£) 




(C.21) 


where a, b and c are positive dimensionless constants (where positive a assumes positive A) 
and an overall factor of 1 /£ 2 comes from the transition to 4D Einstein frame. This implies the 
potential climbs to positive values as i comes in from infinity, eventually reaching a maximum 
and then falling to a minimum before climbing again to infinity as £ —> 0. 
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We can see that the solution and R- describe the maximum of this potential in the low- 
energy theory, and this is why it always occurs at positive values of V e g. It should therefore 
be unstable. It is the solution 1/0.+ and R + that describes the minimum, and whether this 
occurs for positive or negative values of V e g depends on the detailed size of the parameters. 
The condition (4.16) is the condition for the existence of both a local maximum and minimum, 
and when it is not satisfied the potential simply rises monotonically as l falls. (We also see 
that a minimum exists, but always with negative potential, even if A —> 0 provided r A ^ 0.) 


Beyond Rugby Balls 


We here record the properties of the more general bulk solutions appropriate when the source 
branes are not identical. Although the solutions are described constructively here, they may 
also be found by double Wick-rotating the higher-dimensional black-hole solution [30]. 

We start with our standard ansatz, with radial coordinate chosen so that B (x 1/A = 
F(y); that is, 

ds 2 = W 2 (y) g )W dx^dx" + ^ + F(y)c 2 d0 2 . (C.22) 


In these coordinates the (yy) 
becomes 


- (99) Einstein equation, together with Z 

1 ( dTT\ _ 1 fdW\ _ 
B\dp)~AB\Ay)~ ’ 


0 in the bulk, 
(C.23) 


for constant k, and so because AB is y-independent we must have W(y) oc y. This leaves 


ds 2 


—) 9yv &x^dx u + + F(y) cl d6 

VoJ F(y) 


(C.24) 


where yo — or 5 equivalently, k — can be absorbed into a rescaling of the coordinate y. 

The Maxwell field is also given in terms of these metric functions because the Maxwell 
equation is solved by y/—g A ye = Q, which implies 


AyQ — Q Co 



(C.25) 


For definiteness we take the maximally symmetric d-dimensional subspace to be de Sitter (or 
flat) space with Hubble scale H, so R^ = — (d — l)H 2 g^ u and 


n (d) = 


d(d - 1) H 2 

(■ y/vo ) 2 


+ d(d- 1) 



+ - dyF ■ 
y 


(C.26) 


This is to be used in the Einstein equation Tl(d) = —2 k 2 X, with 





(C.27) 
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which implies F(y) must take the form 


F(y) 



= H 2 





2(d—l) 

c 


(C.28) 


where 

_ 2 k 2 A _ k 2 Q 2 

' d(d+ 1)’ d(d-l) 


(C.29) 


and B is an integration constant. 

The constants B and H can be traded for two other parameters, and y_, that define 
the zeroes of F: i.e. F(y±) = 0. Extracting two factors that enforce this vanishing, allows 
one to define 


2(d—l) 


F(y) = A(y + - y)(y - y-)G(y) where G(y) = 1 + ^ 


G n 


n= 1 


with 


G n = < 


given 


Vn , n = 1 

^ yn - ^yo Vn—2 , n = 2,...,d 

{ Vn - ^-yl Vn—2 - f yo +1 Vn-d-1 , n = d + 1, . . . , 2(d - 1) 


yn-=^2 y+y- - * • 

i=0 


77 and 77 are given explicitly in terms of y-t by 


if 2 = 


and B = 


(y+/yo) d+i - (y-/yo ) rf+1 


,(y+/yo ) d_1 - (y-/yo) d_1 _ 

(y+/yo ) 2 - (y-/yo ) 2 
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_(y+/yo) 1_cZ - (y-/yo) 1_d _ 


y± 
. yo 


1 —d / \ 1—d 
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y+ 
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1—d 
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l—d 


c. 


Since H 2 > 0, the above imply 

(y+/yo) 2d - ( y-/yoY 


B > 


(y+/yo) d 1 - (: y-/yo) d 1 


C > 0 V y + > (whenever d > 1). 


(C.30) 

(C-31) 

(C.32) 

(C.33) 

(C.34) 

(C.35) 


As a check, consider the limiting case where y±, yo 1 and H 0. Then, since y n —> n + 1 
in that limit, we find 


C 


d+1 
d - 1 


A , B = 


2d 

d- 1 


A, 


2 , n = 1 

< (n + 1), n = 2,... ,d 

(n + 1) — ^(n — <7), n = d + 1,..., 2(d — 1) 

(C.36) 
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and so 


G(y)->d(d + l)-(d-l ) 2 — 


(C.37) 


Therefore — in this limit — the metric becomes that of a rugby ball with radius, r, and 
defect parameter, a, given by 


r2 = 1 - [(d-Vro]2 ’ “ = “° (1 “ [(d ~ ' 


(C.38) 


where r^ 2 := d(d + l)A = 2k 2 A, which is independent of d. 

From here, we wish to change variables from y to d such that y = y± is identified with 
cosd = ± 1 ; we find that the transformation 


where 


— = W(d ), W(d) := l + r/ cos d 
Vo 


y+ + y- , y+ - y- 

yo = -x- and r) = -■- 

2 y+ + y- 


hlls the bill. Under this coordinate change, we see that 


2 _ y+~y- 


d y = (y+- y)(y - y-) dd = 


sin 2 d d-d 2 


and so the metric becomes 


(C.39) 


(C.40) 


(C.41) 


H 2 = 


(1 ~d) 


ds 2 = W 2 (d) g, lv dx^dx” + r 2 + K(d)a 2 sin 2 d dd 2 \ (C.42) 

= W 2 (d) g^dx^dx" + r 2 {d)(dd 2 + K 2 (d)a 2 0 sin 2 dd9 2 ^) (C.43) 

where we identify Tq 2 := d(d + 1)A = 2/t 2 A and r(d) = ro/K _1 / 2 (d), and where K := d ^ + i) • 
Furthermore, since the expressions for H 2 , B, and y n differ from their unwarped values only 
quadratically in 77 , i.e. 


= W 2 (d) g^v dx^dx 1 ' + r 2 (d) ( dd 2 + K 2 (d)a^ sin 2 d dd 2 


) (l + ^ 3 -^ d 2 ) A - [l + (d — l) V 2 \ C + 0( V 4 ) (C.44) 

(l - ^±^1A+[2 + d(d - l)r f] C + 0(rj 4 ) (C.45) 


y n = (n + l)(l + n< ^ 6 ^ d 2 + C’Cd 4 ) 


(C.46) 


the leading-order corrections to K(d ) arise from the 77 -dependence of the warp factors: 

2(d—1) s 

K(d) ~ 1 - [(d - l)gr 0 ] 2 - f £ 71 G n0 ) g cos d + 0( V 2 ) (C.47) 

' ' n=l ' 


~ 1 — [(d — l)idro] 2 — ( d — - — d[(d — l)i7ro] 2 ) 77 cos-d + 0{rj 2 ) 


(C.48) 
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This means that — if ?? <C 1 — then K± := K($±) is well-approximated by 


K± ^ 1 T ( d - ^ ) rj - [(d - l)Hr 0 ] 2 (l T d V ) + 0{if) 


= (l-[(d-l)Rr 0 ] 2 ) 


fTid-g 


1 - 


^[{d-DHrof 
1 - [(d - l)^r 0 ] 2 


(C.49) 


+ 0(e 2 ). (C.50) 


■=Vh 


Then, since the metric takes the form 

ds± ~ (1 ± 2rj)dsl + ^di? 2 + a^_(d — i?±) 2 d$ 2 ^ 

near each pole (to linear order in e), we find 


r± = 


ro 


K. 


1/2 


iy> 


I”* 


a± = ao-ff± — a 


1 =F ( d - - ) r] H 


+ 0[t?(i? - d±)W] 

+ - d±) 2 ,?7 2 ] 


(C.51) 

(C.52) 

(C.53) 


Rearranging the expression for a± and identifying k 2 L±/(2it) = 1 — a±, we see that 


ot+ + a_ «; 2 T t , r / 9 t- % k 2 /47t 


(L + -L_) + 0[(k 2 L ± ) 2 ]. 

(C.54) 


Furthermore, we see that the warp factor changes by an amount 

K V2ir _ , l 1 - [(d - l|ffr 0 ] 2 


A l C:=M/ W -W(0) = 2 £=((i _ 2/3) 

across the bulk, for H < H cv ; t where 


(T+ - T-) 


1 - 


^[(d-l)Hro]\ 


(C.55) 


[(d l)id crit r 0 ] — 1 3^ ■f - )' -Rcrit — 


(^r) 2b2a ' (a56) 
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